Chiral states in bilayer graphene: magnetic field dependence and gap opening 
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At the interface of electrostatic potential kink profiles one dimensional chiral states are found in 
bilayer graphene (BLG). Such structures can be created by applying an asymmetric potential to the 
upper and the lower layer of BLG. We found that: i) due to the strong confinement by the single 
kink profile the uni-directional states are only weakly affected by a magnetic field, ii) increasing the 
smoothness of the kink potential results in additional bound states which are topologically different 
from those chiral states, and Hi) in the presence of a kink-antikink potential the overlap between 
the oppositely moving chiral states results in the appearance of crossing and anti-crossing points 
in the energy spectrum. This leads to the opening of tunable minigaps in the spectrum of the 
uni-directional topological states. 
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I. INTRODUCTION 

Carbon-based electronic structures have been the fo- 
cus of intense research since the discovery of fullerenes 
and carbon nanotubes 1 . More recently, the production 
of atomic layers of hexagonal carbon (called graphene) 
has renewed that interest, with the observation of strik- 
ing mechanical and electronic properties, as well as 
ultrarelativistic-like phenomena in this condensed mat- 
ter system 2 . In that context, bilayer graphene (BLG), 
which is a system with two Van der Waals coupled sheets 
of graphene, has been shown to have features that make 
it a possible substitute of silicon in microelectronic de- 
vices. The carrier dispersion of pristine BLG is gapless 
and approximately parabolic at two points in the Bril- 
louin zone (K and K'). However, it was found that the 
application of perpendicular electric fields produced by 
external gates deposited on the BLG surface can induce 
a gap in the spectrum, by creating a charge imbalance 
between the two graphene layers 5 ' 6 . The tailoring of the 
gap by an external field may be particularly useful for 
the development of devices 7 ' 8 . It was recently recognized 
that a tunable energy gap in BLG can allow the obser- 
vation of new confined electronic states, which could be 
obtained by applying a spatially varying potential profile 
to create a position-dependent gap analogous to semicon- 
ductor hetero junctions. 

An alternative way to create one dimensional localized 
states in BLG has recently been suggested by Martin et 
al. 9 and relies on the creation of a potential "kink" by 
an asymmetric potential profile (see Fig. 1). Such kink 
potential can also be realized in p-n junctions. They 
showed that localized chiral states arise at the location 
of the kink, with energies inside the energy gap. These 
states correspond to unidirectional motion of electrons 
which are analogous to the edge states in a quantum Hall 
system. From a practical standpoint, the kinks may be 
envisaged as configurable metallic nanowires embedded 
in a semiconductor medium. Moreover, the carrier states 
in this system are expected to be robust with regards to 
scattering and may display Luttinger liquid behavior 10 . 
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FIG. 1: (Color online) (a) Schematic illustration of the non- 
uniformly gated bilayer graphene device for the creation of a 
kink potential. Applied gated voltage to the upper and lower 
layers with opposite sign induces an electric field E e , with 
preferential direction. An external magnetic field B — B z, 
is applied perpendicular to the bilayer graphene sheets, (b) 
Electric field between the two graphene layers, (c) Potential 
on layer 1 (Ui) and layer 2 (L^). 



An additional tool for the manipulation of charged 
states in BLG is the use of magnetic fields. The appli- 
cation of an external magnetic field perpendicular to the 
BLG sheet causes the appearance of Landau levels which 
can be significantly modified by the induced gap, leading 
to effects such as the lifting of valley degeneracy caused 
by the breaking of the inversion symmetry by the elec- 
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trostatic bias 12 ' 13 . Recently the transport properties of 
p-n-p junctions in bilayer graphene were experimentally 
investigated in the presence of a perpendicular magnetic 
field 11 . 

In the present paper we generalize previous work on 
topological confinement in bilayer graphene on three lev- 
els: i) we investigate the effect of smoothing the kink 
potential on the topological states, ii) the effect of a 
perpendicular magnetic field is studied, and Hi) we in- 
vestigate a new system that consists of a coupled kink- 
ant ikink structure. We demonstrate that the latter opens 
a gap in the ID electron states. The paper is organized as 
follows. In Sec. II we present the theoretical formalism. 
The results for a single kink potential profile are discussed 
in Sec. III(A,B). In Sec. IV(A) and Sec. IV(B) we show 
the results for the kink-antink potential, respectively, for 
zero and non-zero magnetic fields. Finally we conclude 
the remarks of the paper in Sec. V. 



II. MODEL 

We employ a two-band continuum model to describe 
the BLG sheet. In this model, the system is described 
by four sublattices, two in the upper (A, B) and two in 
the lower {A' and B') layer 7 . The interlayer coupling is 
given by the hopping parameter t ~ 400 meV between 
sites A and B r . The Hamiltonian around the K valley of 
the first Brillouin zone can be written as 



H = — 





M 2 



(^) 2 




+ 



U{x) 
-U(x) 



(1) 



where tt = v F (p x + ip y ), p XjV = -ihd x , y + eA XiV is the 
momentum operator in the presence of an external mag- 
netic field with A XjV being the components of the vector 
potential A, v F = 10 6 m/s is the Fermi velocity, U(x) 
and —U(x) is the electrostatic potential respectively ap- 
plied to the upper and lower layers. The eigenstates 
of the Hamiltonian Eq. (1) are two-component spinors 
V(x,y) = [^ a (x,y),^b(x,y)] T , where ^a,6 are the enve- 
lope functions associated with the probability amplitudes 
at sublattices A and B' at the respective layers of the 
BLG sheet. Since [H,p y ] = the momentum along the 
^-direction is a conserved quantity and therefore we can 
write: 



ip(x,y) 



— e ik y y 



[if a {x),Lp h {x)} 1 



(2) 



where, k y is the wave vector along the y direction. 

In order to apply a perpendicular magnetic field to the 
bilayer sheet we employ the Landau gauge for the vector 
potential A = (0,Bx,0). The Hamiltonian (1) acts on 
the wave function of Eq. (2) which leads to the following 
coupled second-order differential equations, 

[^7 + (K + P x ')\ 2 Vb = [e~ u{x')]<p a , (3a) 
[JL _ {k > y + $x<)\\ a = [e + u(x')]<p b . (3b) 



where, in the above equations we used dimensionless 
units x' = x/l, k f y = k y l, e = E/t, and u(x') = U(x)/t 
where I = hv F /t = 1.6455 ram, /3 = (eB/h)l 2 (= 0.0041 
for B = 1 T). The step-like kink (see Fig. 1(c)) is mod- 
eled by 



u(x') = lib tanh(a//(5), 



(4) 



where is the maximum value of the gate voltage, in 
dimensionless unit, in each BLG layer. Here, 5 denotes 
the width of the region in which the potential switches its 
sign in each layer. This parameter is determined by the 
distance between the gates used to create the gap. Next, 
we numerically solve Eqs. (3) to obtain the dependence 
of the energy levels on the magnetic field and potential 
parameters. For the case of a sharp kink potential 5 — >• 
and in the absence of a magnetic field , i.e. 5 = 0, Eqs. 
(3) reduces to 

d 2 

[g^j + k' y ] (Pb = [e- u(x')](p a , (5a) 
d 2 

[■faj-K] <Pa = [e + u(x')]<p b . (5b) 

where u(x') = Ub[0(x') — 0(— x')\. We simply decouple 
Eqs. (5) and obtain 



1/2 



(6) 



where, A± = [ — k y 2 ± \/ c 2 — u 2 ] which can be a com- 
plex quantity. The solution for x' < (V ;< ) an d x' > 
( , > ) are given by 



f ±e i\±x' 

X±x' 
-i\±x' 



e 
9±e 



(7a) 
(7b) 



where, f± = (i\± - k' y ) 2 /(e - u h ), and g± = (i\± + 
k y ) 2 /(c + Ub). The above solutions should satisfy the 

asymptotics ^ h ( x ' ~~ ^ 00 ) ® anc ^ ^Pab( x ' ~~ ^ — °°) = 0- 
Matching the solutions and the first derivatives at x' = 
gives a homogeneous set of algebraic equations which in 
matrix form becomes 



0. (8) 



Solutions are found when the determinant of the matrix 
is set to zero from which we obtain the energy spectrum. 
Notice that Eq. (8) leads to four solutions, of which two 
of them, i.e. ±u&, do not satisfy Eqs. (5) and are not 
acceptable. In the limiting case e « M(, we are able to 
obtain an analytical expression for the energy, 




Ub 

a 



[4k' yy fa[u b sin(0/2) + k 2 cos(0/2)] 
: [56k y 8 + 14*4 + ™4ky ~ k 2 e (40k y 4 + 46u b 2 )] 1/2 } . 



(9) 
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FIG. 2: (Color online) (a) Energy levels for a single kink 
profile on bilayer graphene with Ub — 0.25. Dotted curves 
are the numerical results for 5 = 1 and dashed curves are the 
analytical results for 6 = using Eq. (9). The solid red curves 
are the energy levels of a biased BLG. The solid green curve 
(indicated by the symbol (1)) shows a fitted function to the 
numerical results. The lower panels show the wave spinors 
and probability density corresponding to the states that are 
indicated by the arrows (b) and (c) in panel (a). 



where, eo = yk'y +^6 2 , ol = + 7u 2 — 6k y 2 eo and 
6 = t&n -1 (ub/k'y). Solving the above equation for e = 
we find that k' = ±Ju b /V8 (« 0.3 for u h = 0.25). 



Next we consider a sharp kink potential in parallel with 
an antikink potential which are located at x' = — d and 
x' = +d. In this case we have to consider three regions, 
i.e. x' > d (^/(x')), — d < x' < —d (ipjj(x r )) and x' > d 
(ipiu(x f )) and the solutions are given by 

/ i\±x' \ 

M^)±=[ g e ± ^), (10a) 

/ p ±i\±x' \ 
ll>Il{x')± = ( /±eiA± x' J > (10b) 

/ e ~i\±x' \ 

7pIIl{x)±= { 9± e-^') (1 ° C) 

Matching the solutions and their first derivatives at x' = 
±d leads to a set of eight algebraic equations which in 




FIG. 3: (Color online) The probability densities of the topo- 
logical state in k' y = —0.8,-0.6,... which correspond to the 
state indicated by (1) in Fig. (2). The inset shows the prob- 
ability density for k' y — 0.4, 0.6, 0.8 
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where, k,± = ex.p(iX±d) and k± = ex.p(—iX±d). Setting the determinant to zero gives the energy spectrum. 
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FIG. 5: (Color online) Upper panels: applied potential pro- 
file u(x) to the upper (solid blue curve) and the lower layer 
(red dashed curve) for (a) sharp and (b) smooth kink profiles. 
Lower panel: schematic representation of the energy spectrum 
in (E, k y ) space, corresponding to the blue and red dashed po- 
tential profiles near the kink region. Topological states can be 
found in the yellow region. For the smooth profile (b), extra 
bound states can be found in the orange region. 




FIG. 4: (Color online) Upper panel: Energy levels for a single 
kink profile in bilayer graphene with Ub = 0.25 and S = 10. 
The energy states indicated by (a) and (c) are chiral states 
and those indicated by (b),(d),(e) and (f) are the extra-bound 
states. Lower panels: Real parts of the wave spinors and the 
corresponding probability density for the two first electron 
and hole energy levels at k' y — as indicated in the upper 
panel. 



FIG. 6: (Color online) The carrier velocity in the single kink 
profile for the energy levels which are indicated by (1), (2), ... 
in Fig. 4. The gray solid curves are the fitted functions to 
curve (5) and (6). 



III. SINGLE KINK 
A. Influence of the smoothness of the kink profile 



In the general case of 5 ^ we solve the set of second 
order differential Eqs. 3(a,b) numerically, using the finite 
difference technique. Figure 2(a) shows the spectrum for 
a single potential kink as function of the wavevector along 
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FIG. 7: (Color online) Number of additional bound states as 
function of the width of the interface 5 for Ub = 0.1, 0.25 and 
0.5. 
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FIG. 8: (Color online) (a,b) Contour plot of the transmission 
T (in logarithmic scale) through a single kink structure with 
L x 20 (« 32 nm), u h 0.25 for (a) 8 0.1 and (b) 5 = 4. 
(c) Conductance G/Go vs Fermi energy of a single kink profile 
for the same parameters as (a,b). 



the kink for zero magnetic field. We consider a relatively 
sharp kink, i.e. 5=1, and compare the numerical re- 
sults with the analytical solution (dashed black curves) 
form Eq. (9) for the case of a sharp profile (S = 0). 
The shaded region corresponds to the continuum of free 
states. The solid red curves correspond to the energy 
levels of a biased BLG which can be obtained using Eq. 
(1) as, 



e = ± 



The dotted horizontal lines correspond to e = ±ui> = 
±0.25 and e = 0. These results are valid in the vicin- 
ity of a single valley (K) and show that the topo- 
logical states have a unidirectional character of prop- 
agation, i.e. they are chiral states 9 , with positive 
group velocity. The topological levels can be fitted 

to e = a^j (k' y — ko) A + (u^/a) 2 — 1/2 with a = 0.5 and 

fco = —0.65 being the fitting parameters (see green solid 
curve) . For localized states around the K' valley, we have 
Ex f (k y ) = —Exiky) and the charge carriers move in the 
opposite direction. In order to consider the energy levels 
for the K' valley, u{x) in Eqs. (3) should be replaced 
with —u(x). Then using the transformations e — » — e and 
(p a —> —(p a (or tpi — >> —ifb) leads to the same equation as 
for the K valley. Thus the E K (k) = —E' K (k) symmetry 
remains even in the presence of an uniform perpendicular 
magnetic field (i.e. ^ 0). Notice that the wavespinors 
corresponding to the K and the K' valleys are related to 
each other by (p^ — > —(fa ( or ^ ~~ ^ 
sign of the other component does not change. 

Panels (b) and (c) of Fig. 2 present the real parts of 
the spinor components and the probability density for 
the states indicated by the arrows (b) and (a) in panel 
(a), corresponding to k! y = —0.28 (b) and k' y = 0.28 (c). 
These electron states are localized at the position of the 
potential kink. Notice that the solutions of Eqs. (3) are 
related by the transformations ip a — > — <£b, (fb —> Pa, 
k y — >• — k' y and e — » — e and consequently for e = the 
solutions in Figs. 2(b) and 2(c) have the same probability 
distribution. For the case of k' y = the solutions of 
Eq. (9) are e± = ±Ub/V2 which result in the following 
wavespinors, 



(fif' ) while the 



^ = (lT^)[e-^'^Te 



(1 =F V^) r A 2 e -i(A + x' + ^) A 2 e i 

(e± - u b ) + 
-1 



(13a) 
(13b) 

(A -*' + *>],(13c) 



e± + u b 



[ x 2 +e i(X +X ' + V T> 2_ e -iX_ X 'y 



(13d) 



where, X± = 7(1 ± i) with 7 being y/u^/S. Notice that in 
the above equations e ±%x±x leads to an oscillating con- 
tribution e ±^( A ±)^ with an evanescent e^^ x± ^ x part. 
The oscillating part is strongly damped and therefore 
Eqs. (13) corresponds to localized wavespinors. Expand- 
ing Eqs. (13) around x' = we obtain for the second 
derivative of the wavespinors, 



d 2 

dx' 2 
d 2 
dx' 2 



St[<p<(x' 



0)] = (1± V2) (2 + ^2)7, (14a) 

o)] = ii^Mh! < o. (14b) 



e± - u h 



K 4 



(12) 



This indicates that 5£[</? a ](5£[^>]) has its maximum value 
located at x' < 0(x' > 0) for e + = u^/V^ while the 
opposite is found for e_ = —u^/V^ which is also evident 
from Figs. 2(b,c). 
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In Fig. 3 we show the probability densities correspond- 
ing to one of the topological branches (for the state which 
is labeled by (1) in Fig. 2(a)) at several k' y values. As 
shown in the inset of Fig. 3 for those k' y values where the 
topological state merged with the continuum spectrum 
the carries are no longer confined by the kink potential. 

Next we increase the smoothness of the kink potential 
and investigate how the energy spectrum changes. In 
Fig. 4 the energy levels as function of k y are shown for 
the smooth kink profile 5 = 10, where in addition to the 
chiral states several branches are seen which are split off 
from the continuum. In order to understand the physical 
origin of those new states we show in the lower panels 
of Fig. 5 a cartoon of the low energy spectra for the (a) 
sharp and (b) smooth profiles where, the chiral states ap- 
pear in the yellow regions and those additional states are 
found in the orange region. Increasing the smoothness of 
the kink potential leads to the creation of a region below 
the energy gap which allows for carriers to be confined 
near the kink. Therefore, extra bound states can be cre- 
ated in the orange region (lower panel in Fig. 5(b)). The 
wavefunctions for k y = of the two chiral states and the 
new bound states are shown in the lower panels of Fig. 4. 
The new bound states are also bound in the x-direction 
near x = but the electron states are more extended and 
have a clear nodal character near x = 0. 

Figure 6 shows the velocity of the carriers for the states 
which are indicated by (1),(2),... in Fig. 4. The chi- 
ral states ((5), (6)) are only shown for the K valley and 
they have positive velocity. The curves (5), (6) can be 

fitted to v/vp ~ (ky — k' ) 2 + b 2 (see the solid gray 

curves) with a = 1.8,6 = 0.15 being the fitting pa- 
rameters and k f = ±0.08 corresponds to the minimum 
point in the curves (5) and (6). Notice that the extra 
bound states ((2), (3), (4)) have a slightly nonzero ve- 
locity at k' y = which is a consequence of the asym- 
metric energy dispersion as seen in Fig. 4. Curve (1) 
corresponds to the energy spectrum of a biased BLG 
which is given by Eq. (12) and results in the veloc- 
ity v/v F = [de/dk'y 



u 2 which is zero for 



k' y = in a biased BLG (black solid curve in Fig. 6). 

As mentioned before for smooth kink potentials addi- 
tional ID bound states appear and the number of these 
bound states can be related to the height of the gate volt- 
age lib and the smoothness(^) at the interface. Figure 7 
shows the number of these extra bound states for three 
different values as function of the width 8. The first 
bound state for = 0.1,0.25,0.5 appears respectively 
at 5 « 6, 4, 2 in the absence of magnetic field. Notice 
also that for fixed 5 the number of extra bound states in- 
creases with Ub in agreement with the qualitative picture 
shown in Fig. 5(b). 

We also calculate the transmission of an electron 
through the kink structure in a system of size —L x /2 < 
x < L x /2 and —L y /2 < y < L y /2. No bias nor magnetic 
field is assumed in the x < —L x /2 and x > L x /2 regions. 
We assume that L y » L x and the electrons are free to 
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FIG. 9: (Color online) Energy levels of a single potential kink 
profile as function of the external magnetic field with Ub = 
0.25 and S = 10 for (a) k y l = and (b) k y l = 0.15. The 
full blue curves are the topological states and the dashed red 
curves are the extra bound states. 



move in the ^/-direction whereas, they are confined in the 
x-direction. Associated with each real A± there are two 



right(left) propagating modes, tp±(ij)± 
by Eqs. (7). In the region I (x < 
dent right-traveling modes ip± can be reflected into two 
left-traveling modes 



which are given 
-L x /2) two inci- 



(15) 



where, t± (r±) are the transmission (reflection) ampli- 
tudes. The propagating modes in region I can also 
be transmitted to region III (x > L x /2) in the right- 
traveling modes, 



(16) 



The wevefunctions in regions I and III can be connected 
by the transfer matrix M where at the kink-potential 
boundaries we have 



* I ±(-L x /2)=M^ I (L x /2). 



(17) 



The transmission (or reflection) amplitude can be found 
by substituting Eqs. (15) and (16) in the above equation. 
The four transmission amplitudes t± for given e and k y 
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FIG. 10: (Color online) (a) Energy spectrum of a single kink 
profile in bilayer graphene as function of the cyclotron orbit 
coordinate X c for B = 7 T, Ub = 0.25 and S = 1. The dashed 
gray curve shows a fitted function to the numerical results. 
The solid black and green curves are respectively the potential 
in the upper and lower layer. Solid red lines corresponds to 
the Landau levels of a biased BLG. The solid purple curves 
show fitted functions (given in the text) to the position of the 
resonances. (b,c) The probability densities for the indicated 
points by red full circles in the spectrum. 



FIG. 12: (Color online) ^/-component of the persistent cur- 
rent corresponding to the topological state of a single kink 
potential as function of the x direction for zero magnetic field 
(dashed curves) and B — 5 T (solid curves) with k' y — and 
Ub = 0.25. Blue curves display the current density for 5 = 1 
and black curves are the corresponding results for 5 = 10. 
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FIG. 11: (Color online) The same as Fig. 10 but for S = 10. 
The solid black and green curves are the potential, respec- 
tively, in the upper and lower layer. 



FIG. 13: (Color online) Density of states (DOS) of a single 
kink potential for S = 1 (solid and dashed curves) and 5 = 10 
(dotted and dashed-dotted curves) with Ub = 0.25 for two 
different magnetic field values. 

can be combined in the transmission matrix 

t(e,k y )= (*t (18) 

The total transmission amplitude T is given by 17 T = 
Tr(ttf). The two-terminal conductance of such an asym- 
metric potential profile in bilayer graphene can be calcu- 
lated using the Landauer formula which is given by 18,19 

G = G f T(E F ,k' y )dk' y (19) 

Here, Go = (2e 2 /h){L y \EF\/iihvF) is the conductance 
unit per valley and per spin. In Figs. 8(a,b) we show 
a contour plot of the transmission probability in loga- 
rithmic scale for the kink structure with L x = 20 (in di- 
mensionless unit). The transmission probability has the 
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FIG. 14: (Color online) Conductivity a yy as function of the 
Fermi energy Ef for a single kink potential with Ub — 0.25. 
Other parameters are indicated in the figure. 

symmetry T(k y ,E) = T(—k yi —E). The conductance as 
function of the Fermi energy for the single kink profile is 
shown in panel (c) for 6 = 0.1 (blue solid curve) and 6 = 4 
(red dashed curve). For the case 6 = 4 the smoothness 
of the potential at x = leads to a higher transmittance 
and consequently a higher conductance around e ~ u b 
(see panel (b) and the dashed curve in panel (c)). 



B. Magnetic field dependence 

Dependence of the energies of the ID bound states on 
an external magnetic field is shown in Fig. 9 for (a) 
k f y — and (b) k y = 0.15. In order to show the ef- 
fect of a magnetic field on the chiral states (blue solid 
curves) and the other localized bound states (red dashed 
curves) we present the results for a smooth potential (i.e. 
6 = 10). It is seen that the chiral states are very weakly 
influenced by the magnetic field. This is a consequence 
of the strong confinement of these states in the kink po- 
tential (see Fig. 3 and Fig. 4(a,c)). In a semiclassical 
view, the movement of the carriers is constrained by the 
kink potential and that together with the unidirectional 
propagation, prevents the formation of cyclotron orbits. 
For the energy levels above the chiral states, the energy 
values increase as the magnetic field increases, because 
of the weaker confinement of these states as is apparent 
from Figs. 4(b,d,e,f). 

Figure 10 shows the spectrum of a sharp (6 = 1) sin- 
gle kink potential in the presence of an external mag- 
netic field B = 7 T as function of the orbit center 
X c = k y l 2 B /l where Ib — ^JhjeB is the magnetic length. 
The solid lines represents the applied kink potential to 
upper (black) and lower (green) layer. The results show 
that the topological states are practically not affected by 
the magnetic field. The free energy region (i.e. |e| > u b ) 
in the absence of magnetic field now is replaced with Lan- 



dau levels (the solid red lines are the Landau levels of a 
biased bilayer graphene). In some region the Landau lev- 
els are influenced by the kink potential and anti-crossings 
appear in the low energy spectrum. Some of these anti- 
crossings are situated along the extension of the topolog- 
ical states into the |e| > region. In addition of these 
anti-crossings the Landau levels display some resonances 
along the energy levels of a biased BLG (red solid curves 
in Fig. 2(a)) which can be linked to the edge effects of the 
potential profile. The position of the resonances can be 
fitted to e = ay/(X c - X ) 4 + (e n /a) 2 where a = 0.0006 
and Xq = 9.25 are fitting parameters and e n indicates 
the nth Landau level of a biased BLG (see solid pur- 
ple curves). Also the topological levels can be fitted to 
e = (X c — X ) 4 + (ub/a) 2 (dashed gray curve) with 
a = 0.0003 and Xo = —26. Panels (b,c) show the proba- 
bility densities for the points that are indicated by full red 
circles in the energy spectrum. For the points on the pur- 
ple solid curves (2a,3a) the distribution of the carriers by 
the magnetic field is influenced by the small confinement 
by the interface potential (see solid and dashed curves in 
panel b) . The probability density for the point on the fit- 
ted curve along the topological level (2b) shows a higher 
peak at the kink interface (x = 0) indicating that the 
kink potential acts as an attractive potential (solid curve 
in panel (c)). The other probabilities are clearly those of 
free electron LL. The result for a smooth kink potential 
5 = 10 is shown in Fig. 11 where the energy value of the 
extra bound states are increased by the magnetic field 
and the topological levels are practically not affected by 
the magnetic field. 

The localization of the states is reflected in the position 
dependence of the current. The current in the ^/-direction 
is obtained using 

j y = iv F [9\d x a y -dya x )9 + 9 T (d x ay + d y a x )^ ( 20 ) 

where \I/ = [ip a (x,y) , , y)) T '• By substituting 

ty(x,y) = e lkyV [(p a (x), (pb(x)] T we have 

j y = 2v F [Re{(p* a d x (p b - (ftd x (p a } + 2k y Re{<pZ<p b }] . (21) 

The x-component of the current vanishes for the confined 
states. In Fig. 12, the ^/-component of the persistent cur- 
rent for a sharp (blue curves) and smooth (black curves) 
single potential kink profile is shown as function of the x 
direction without magnetic field (solid curves) and in the 
presence of the magnetic field (dashed curves) . In the ab- 
sence of a magnetic field the current is localized around 
x = for both sharp (6 = 1) and smooth (6 = 10) po- 
tentials. For a smooth profile the wavefunction of the 
topological states and consequently also the current den- 
sity profile is broadened (compare Figs. 2(b,c) with Figs. 
4(a,c)). A magnetic field shifts the density profile slightly 
to the right (see the inset of Fig. 12) due to the Lorentz 
force and there is also a very small narrowing of the cur- 
rent distribution. 

Next we consider the density of states (DOS) for the 
kink potential. The number of k-states per unit energy 
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FIG. 15: (Color online) Left panels: Energy levels of a kink-antikink profile in bilayer graphene with Ub — 0.25, 8 = 1, d = 15 
for (a) B = T and (b) B = 3 T. The black dashed curves display analytical results as obtained from Eq. (11). Right panels: 
Real part of the wavespinors and the corresponding probability density for the points that are indicated in the panels (a,b). 



is given by 

D{E) = § E/ d Mfc " (22) 

where Do = (7wf) _1 . To calculate the DOS numerically 
we introduce a Gaussian broadening, 

,5(e-e n)fe; ) exp[-^^-], (23) 

where T is the broadening which is taken as Y = 0.005 in 
our calculations. Figure 13 shows the DOS as function 
of Fermi energy Ep in the absence and presence of an 
external magnetic field for sharp (5 = 1) and smooth (5 = 
10) kink potentials. For a sharp profile the topological 
levels contribute an almost constant value to the DOS for 
|e| < Ub even in the presence of an external magnetic field. 
For the smooth profile, peaks corresponding to the non- 
topological levels appear in the DOS and note that only 
these peaks are shifted in the presence of a magnetic field 
while the DOS of the topological states are not affected 
by the magnetic field. 

We now turn to the transport properties of a kink po- 
tential and look at the influence of the topological states 
on the conductivity in the ^-direction (<J yy ). For elastic 



scattering the diffusive conductivity a yy is given by 20 , 
ff » = £ / dk >ljn, K {l ~ /»,*;)• (24) 

Here T is the temperature, v Uyy = de n /dk' y is the elec- 
tron velocity, / n>fe = l/[exp ((E F - E n # y )/K B T) + 1] is 
the equilibrium Fermi-Dirac distribution function, and 
r is the momentum relaxation time. For low tempera- 
tures we assume that r is approximately constant, eval- 
uated at the Fermi level (r « rp), and replace the prod- 
uct f n ,k' y 0- — fn,k' y )/kBT by the delta function given in 
Eq. (23). The results are presented as function of Ep in 
Fig. 14 in the units of <Jo = ^vpTp jTiL y for both sharp 
(6 = 1) and smooth (S = 10) potentials with B = T 
and B = 5 T. Due to the robust confinement of the 
topological levels the conductivity is constant in the en- 
ergy gap even for a non-zero magnetic field (solid blue 
curve for B = T and red dashed curve for B = 5 T). 
The extra localized levels in the case of 5 = 10 lead to 
an increasing conductivity as function of Ep. Note that 
in the presence of an external magnetic field some of the 
additional electron(hole) states are shifted up(down) in 
energy (see Figs. 4 and 11) which results in smaller a yy at 
the |e_p| < Ub region compared to the conductivity in the 
absence of magnetic field (black dotted-dashed curve). 
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IV. KINK-ANTIKINK 

A. Zero magnetic field 

Next we consider a potential profile with a pair kink- 
antikink. The kink-antikink potential is modeled by, 

u(x') = u b [tcmh( X ~ d ) - tanh( X ~^ d ) + l] (25) 

where, d is the distance between the kink and the an- 
tikink in units of I. The spectrum of the localized states 
in the absence of a magnetic field is shown in Fig. 15(a) 
for lib — 0.25, S = 1 and d = 15. The black dashed curves 
are the analytical results for 5 = which are obtained us- 
ing Eqs. (11). Note that there are only two chiral states 
per kink which leads to the appearance of crossing points 
in the energy spectrum (at E = and k' y — 0). The 
spinor components and probability densities associated 
with the points indicated inside the circle in Figs. 15(a) 
are shown in the panels (la, 2a,. ..,5a). In the absence 
of a magnetic field and for the points around the energy 
level crossing the carriers are strongly confined at either 
the position of the kink or antikink. The wavefunction 
corresponding to an energy at the crossing point (panel 
5a) is localized at both the kink and antikink. 

Next we investigate smooth potential kink profiles. 
In Figure 16(a) the energy spectrum of a smooth kink- 
antikink profile (i.e. 5 = 10) is presented for zero mag- 
netic field. As in the case of the single kink profile addi- 



tional bound states appear in the energy spectrum. The 
overlap between theses states leads to the appearance of 
crossing points in the energy spectrum. The wavespinors 
and the corresponding probability density for the points 
indicated by arrows in panel (a) are shown in panels 
(la, 2a, 3a, 4a). In the absence of a magnetic field and for 
k' y = the states are localized at both kink and antikink 
(panels la, 2a and 4a) whereas, panel 3a shows that the 
confinement tends to the kink or antikink at k' y ^ 0. 

Decreasing the distance between the kink and antikink 
generate an unperfect kink-antikink profile 21 . This pro- 
file is illustrated in Fig. 17(a). The energy spectrum of 
such a profile is shown in Fig. 17(b) for B = T, 5 = 1 
and d = 3(« 5 nm). The analytical results (obtained 
from Eq. (11)) for S = are shown by the black dashed 
curves. Now the crossing points in the energy spectrum 
for the case of d = 15 (see Fig. 15(a)) are replaced with 
anticrossings and an energy gap E g appears in the en- 
ergy spectrum. The positions of these minigaps move 
when we increase the magnetic field as in apparent from 
Fig. 17(c). The panels (lb, 2b, 3b, 4b) show the real parts 
of the wavespinors and corresponding probability density 
for the indicated points in Fig. 17(b) by red arrows. Note 
that due to the decreasing distance of the kink and an- 
tikink the carriers can be localized between the kink and 
antikink. 

Figure 18(a) displays the energy spectrum of a smooth 
(S = 10) kink-antikink potential with d = 3 for B = T. 
Now the kink and antikink are close to each other and 



11 



(a) 







\ / 


\ / 
\ / 

\[ 




-6 A 


* 6 

y \ 

' v» 


► 

x/l 







/ y 




<Pb'- 


U 1b 



-0.6 -0.4 -0.2 0.2 0.4 



ky I 



/ \ 

/ 




/ 










x/l 



10 10 

x/l 



FIG. 17: (Color online) Energy levels of a sharp kink-antikink profile with Ub — 0.25, d — 3(~ 5 nm) and 5 — 1 (the profile is 
depicted in (a)) for (b) B = T and (c) B = 3 T. The black dashed curves in panel (b) display analytical results as obtained 
from Eq. (11). Right panels show the real part of the wavespinors and the corresponding probability density for the points at 
the energy spectrum that are indicated by arrows in panels (b,c). 



the smoothness of the potential leads to extra localized 
levels. Therefore the crossing and anticrossing points be- 
tween the additional bound states are seen to disappear 
and the energy gap between the topological levels is in- 
creased. The magnitude of the energy gap E g depends on 
the width of the interface region (5, the maximum value of 
the potential and the distance between the kink and 
antikink. This is shown in Fig. 19, where E g is plot- 
ted as function of 1/5, 5 and d respectively in panels (a), 
(b) and (c) in the absence of magnetic field (blue solid 
curves). As shown in panels (a,b) the energy gap is an 
increasing function of and S. When S increases the 
first energy level at the spectrum changes from a Mex- 
ican hat shape to a parabola. Therefore, E g increases 
with increasing S (compare the potentials illustrated in 
Figs. 17(a) and 18(a)). Increasing the distance of the 
kink and antikink results in perfect unidirectional states 



and the gap disappears (panel 19(c)). 



Next we consider the transmittance of a kink-antikink 
potential. In Fig. 20 we show a contour plot of the trans- 
mission probability (in logarithmic scale) for the kink- 
antikink structure with L x = 24 for (a) 6 = 0.1 (sharp) 
and (b) S = 4 (smooth) potentials. The results show a 
nonzero region for the transmittance below the gap where 
the topological levels corresponding to the kink and an- 
tikink cross each other. The conductance as function of 
Fermi energy is plotted in Fig. 20(b). A small region of 
transmittance appears in the energy gap due to the chi- 
ral states that appears as small peaks in the conductance 
(see the inset of panel (c)). 
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FIG. 18: (Color online) Energy levels of a smooth kink- 
ant ikink profile on bilayer graphene with Ub = 0.25, d = 3 
and 5=10 for (a) B = T and (b) B = 3 T. 




FIG. 19: (Color online) The energy gap (indicated in Fig. 
17(b)) as function of (a) maximum amplitude of the gate volt- 
age Ub, (b) the width of kink profile S and (c) the position of 
the kink and antikink d. Other parameters are shown in the 
figures. 
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FIG. 20: (Color online) (a,b) Contour plot of the transmis- 
sion T (in logarithmic scale) through a kink-antikink structure 
with the length L = 24 (« 40 nm) for Ub — 0.25 and d = 6 
with (a) S = 0.1 and (b) S = 4. (c) Conductance G/Go vs 
Fermi energy of a kink-antikink potential for the same param- 
eters as (a,b). 
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FIG. 21: (Color online) Energy spectrum of a kink-antikink 
profile in bilayer graphene as function of the cyclotron orbit 
center X c with Ub — 0.25, d = 15, B = 7 T and 5 = 1. 
Solid red lines correspond to the Landau levels of a biased 
BLG. The gray solid curves show fitted functions (given in 
the text) to the numerical results. The solid black and green 
curves describe the potential, respectively, in the upper and 
lower layer. 
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FIG. 22: (Color online) Energy levels of a sharp kink-antikink 
profile (5 = 1) as function of external magnetic field for Ub — 
0.25 with (a) k' y = and (b) k' y = 0.2. 
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FIG. 23: (Color online) y-component of the persistent current 
of a kink-antikink profile in bilayer graphene as function of the 
x-direction for E — Ef and for values of the magnetic field 
that are indicated by (1), (2), ... in Figs. 22 (a, b). 



B. Magnetic field dependence 

Figure 15(b) shows the kink-antikink energy levels in 
the presence of an external magnetic field (B = 3 T). The 
results show a shift of the four intra-gap energy branches 
as the magnetic field increases. In addition, the con- 
tinuum of free states at zero magnetic field (shadowed 
region in Fig. 15(a) is replaced by a set of Landau levels. 
The spinor components and probability densities associ- 
ated with the points indicated inside the circle in Figs. 
15(b) are shown in the panels (lb, 2b,... 5b). For non-zero 




FIG. 24: (Color online) Density of states (DOS) of kink- 
antikink potentials for (a) S = 1 and (b) 6 — 10 with 
Ub = 0.25. The results are presented for d = 3 and d — 15 
with B = T and B = 5 T. 



magnetic field the states show a shift of the probability 
density towards the region between the kink and the an- 
tikink. This is caused by the additional confinement due 
to the magnetic field. 

The energy levels of a smooth kink-antikink profile (i.e. 
S = 10) in the presence of a perpendicular magnetic field 
is presented in Fig. 16(b). Now the crossings points (in 
the case of B = T) changed into anticrossings. In the 
inset of Fig. 16(b) an anti-crossing is enlarged. Due 
to the strong confinement of the potential the magnetic 
field can only lead to a shift up in energy of the localized 
chiral states. The wavespinors and the corresponding 
probability density for the points indicated by arrows in 
panel (b) are shown in the panels (lb, 2b, 3b, 4b). In the 
presence of an external magnetic field and at the cross- 
ing points of the topological states (panels lb, 2b) due 
to the strong confinement by the potential the magnetic 
field can only affect weakly the electrons. At the first 
anticrossing (panels 3b and 4b) which arises form the 
overlap of the first bound states in the kink and antikink 
potentials, the electrons are confined closer to the center 
of the potential. 

The energy levels for a sharp (S = 1) kink- antikink 
potential with d = 3 is presented in Fig. 17(c). The 
crossings which appeared in the energy spectrum due to 
the overlap of the extra bound states in the absence of 
magnetic field (see Fig. 17(b)) now are replaced with 
anti-crossings and the energy gap E g between the kink 
and antikink states is shifted up in energy due to the con- 
finement by the magnetic field. Panels (lc, 2c, 3c, 4c) show 
the wavespinors and probability density for the points in- 
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dicated by arrows in panel (c). The energy spectrum of 
a smooth kink-antikink potential with d = 3, 5 = 10 and 
in the presence of an external magnetic field B = 3 T 
is shown in Fig. 18(b). In the presence of a magnetic 
field the energy gap is shifted and the symmetry of the 
spectrum around E = for B = T is broken (see panel 
(a)). The energy gap in the presence of a magnetic field 
(B = 5 T) is shown in Fig. 19 as red dashed curves. 
Notice that an external magnetic field only shifts up the 
energy gap and the gap size remains constant. 

The energy spectrum of a kink-antikink potential is 
shown in Fig. 21 as function of orbit center X c for S = 1, 
Ub — 0.25 and d = 15. The kink-antikink potential is de- 
picted in the figure by the dashed curves. Such as for the 
single kink potential the topological levels can be fitted 
to e± ~ a^/ (X c db Xq) 4 + (lib/ a) 2 (see gray solid curves) 
where — (+) corresponds to the kink(antikink) branches 
(a — 0.003 and X = 41 are the fitting parameters). 
Now the Landau levels above the gap are affected by the 
kink-antikink potential where anti-crossing points appear 
along the topological levels. The solid red lines are the 
Landau levels in a biased BLG. 

Figure 22 shows the dependence of the energies on the 
external magnetic field for (a) k' y = and (b) k' y = 0.2. 
The branches that appear for |e| > 0.25 correspond to 
Landau levels that arise from the continuum of free 
states. For the kink-antikink case, however, the over- 
lap between the states associated with each confinement 
region allows the formation of Landau orbits. Therefore, 
the proximity of an antikink induces a strong dependence 
of the states on the external field. 

Figure 23 shows plots of the ^-component of the cur- 
rent density as function of x for the states labeled (1) 
to (6) in panels (a) and (b) of Fig. 22. It should be 
noticed that a non-zero current can be found for E = 
and k' y ^ 0, as can be deduced from the dispersion re- 
lations. For k' y = the results presented in Fig. 23(a) 
show a persistent current carried by electrons localized 
at each kink region, irrespective of the direction of as 
exemplified by the states (1) and (2) which correspond 
to opposite directions of the magnetic field. For non-zero 
wave vector, however, as shown in panels (b) and (c), the 
current is strongly localized around one of the potential 
kinks. In Fig. 23(b), the current density curve shows an 
additional smaller peak caused by the strong magnetic 
field (B & 7.5 T) where, the carriers can also be confined 
closer to the center. 

The density of states of the topological states for 
(a) 5 = 1 and (b) 5 = 10 kink-antikink potential is 
shown in Fig. 24 with d = 3 and d = 15. The results 
show additional peaks for a sharp kink-antikink with 
d = 3 which is due to the splitting of the topological 
levels. Note that the energy gap leads to a zero density 
at Ep = for zero magnetic field (blue circles in (a)) 



while shifting the gap in the presence of a magnetic field 
results in a non-zero DOS at Ep = (red diamonds 
in (a)). For the smooth profiles the non-topological 
ID states lead to the appearance of additional peaks 
in the DOS (panel (b)) that shift with the magnetic field. 



V. CONCLUDING REMARKS 

In summary we obtained the energy spectrum, the den- 
sity of states, the transmission and conductivity for carri- 
ers moving in BLG in the presence of asymmetric poten- 
tials (i.e. kink and kink-antikink profiles) in each layer 
of the BLG. Uni-directional chiral states are localized at 
the location of the kink (or antikink). By controlling 
the gate voltages and/or the smoothness of the kink pro- 
file the number of one-dimensional metallic channels and 
their subsequent magnetic response can be configured. 

The effect of an external magnetic field perpendicular 
to the bilayer sheet was investigated. We found that the 
influence of the magnetic field is very different for single 
and double kinks. Due to the strong confinement by the 
kink potential, the topological states are weakly affected 
by the magnetic field in the case of a single kink profile. 

Changing the sign of the kink potential smoothly (i.e. 
broadening the kink potential) leads to extra bound 
states which have a very different behavior as compared 
to the uni-directional topological states. First, these 
states are no longer uni-directional and they have a quasi- 
1D free electron- type of spectrum which is asymmetric 
around k y = 0. Second, they are less strongly localized 
at the kink of the potential as compared to the chiral 
states and their probability distribution appears as those 
of excited states of the chiral state. 

In the case of parallel kink-antikink profiles apparent 
crossings of the energy levels are found in the spectrum. 
Decreasing the distance between (and/or smoothing) the 
kink-antikink profiles turns into anti-crossings. It opens 
a gap in the topological state spectrum. This allows for 
a robust ID system having a tunable minigap. 
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